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Pd43NiioCu27P20 has been investigated in its equilibrium liquid state with incoherent, inelastic 
neutron scattering. As compared to simple liquids, liquid PdNiCuP is characterized by a dense 
packing with a packing fraction above 0.5. The intermediate scattering function exhibits a fast re- 
laxation process that precedes structural relaxation. Structural relaxation obeys a time-temperature 
superposition that extends over a temperature range of 540 K. The mode-coupling theory of the liq- 
uid to glass transition (MCT) gives a consistent description of the dynamics which governs the 
mass transport in liquid PdNiCuP alloys. MCT scaling laws extrapolate to a critical temperature 
T c at about 20 % below the liquidus temperature. Diffusivities derived from the mean relaxation 
times compare well with Co diffusivities from recent tracer diffusion measurements and diffsuivities 
calculated from viscosity via the Stokes-Einstein relation. In contrast to simple metallic liquids, 
the atomic transport in dense, liquid Pd43NiioCu27P20 is characterized by a drastical slowing down 
of dynamics on cooling, a q~ 2 dependence of the mean relaxation times at intermediate q and a 
vanishing isotope effect as a result of a highly collective transport mechanism. At temperatures 
as high as 2 x T c diffusion in liquid Pd43NiioCu27P20 is as fast as in simple liquids at the melting 
point. However, the difference in the underlying atomic transport mechanism indicates that the 
diffusion mechanism in liquids is not controlled by the value of the diffusivity but rather by that of 
the packing fraction. 

PACS numbers: 61.25.Mv,61.20.Lc,61.12.-q 



I. INTRODUCTION 

We investigate microscopic dynamics in liquid 
PdNiCuP melts with inelastic neutron scattering. Our 
results show that PdNiCuP can in good approximation 
be regarded as multicomponent hard-sphere like system. 
As compared to liquid alkali-metals the packing in liquid 
PdNiCuP is much more dense. As a consequence, dynam- 
ics in liquid PdNiCuP can not be described by concepts 
developed for simple liquids. Instead, atomic transport in 
PdNiCuP melts is in excellent accordance with concepts 
developed in the context of glass formation. 

Because of the short range nature of the interatomic 
potential with a strong repulsive core, the potential in 
simple metals can in good approximation be condensed 
to an effective hard sphere radius R. Therefore, the pack- 
ing fraction ip= |tt n R 3 (n is the number density in unit 
volume) is an important parameter for the discussion of 
the transport mechanism JjJ]. Due to the low compress- 
ibility of liquid metals the hard sphere radius is fairly 
temperature independent Liquid alkali metals are a 
paradigm of hard-sphere like fluids and their microscopic 
dynamics have been investigated in great detail || . At 
low packing fraction, e.g. at temperatures well above the 
melting temperature, atomic transport is dominated by 
binary collisions. This reflects itself in a Gaussian line 
shape of the quasielastic neutron scattering signal in the 
free particle limit towards large wavenumbers q. With in- 
creasing packing fraction, i.e. by approaching the melting 
temperature T m , fluid dynamics play a more important 
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role and towards small but finite q the quasielastic line 
is well approximated by a Lorentzian function || ^| . Al- 
kali melts exhibit a packing fraction of about 0.4 at the 
melting point. 

Colloidal suspensions B and molecular dynamics sim- 
ulated binary metallic melts ]?], |), that have a signif- 
icantly higher packing fraction above 0.5, exhibit mi- 
croscopic dynamics that is not governed by binary colli- 
sions and is typical for glass forming, molecular liquids 
H |llj . Their dynamics has been described by predic- 
tions of the mode coupling theory (MCT) of the liquid 
to glass transition fll^ : Atomic transport slows down by 
orders of magnitude on a small change in temperature or 
pressure. The slowing down of the dynamics goes along 
with a spread of the quasielastic line as compared to the 
Lorentzian found in simple liquids corresponding to a 
stretching in time of correlation functions over a wider 
time range than expected for an exponential decay. 

Mode coupling theory gives a microscopic explanation 
for this behaviour in terms of a non linear coupling of den- 
sity fluctuations caused by feedback effects in the dense 
liquid. Atomic transport is envisioned as a highly collec- 
tive process. At a critical packing fraction ip c , or a critical 
temperature T c respectively, MCT predicts a change in 
the transport mechanism from liquid-like flow to glass 
like activated hopping processes. MCT calculations for a 
hard sphere system exhibit a critical packing fraction of 
0.525 [13]. In a colloidal suspension <p c ~ 0.58 has been 
found pj. Whereas in alkali melts the packing fractions 
are well below 0.5, new multicomponent alloys on Zr- |l4|] 
and Pd-basis jl5|, [l(| [l?]] at quasieutectic compositions 
exhibit packing fractions at the their liquidus tempera- 
ture Tu q of some 0.52 jl8|, [l9|. These alloys known for 
their bulk glass forming ability, exhibit viscosities at their 
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melting temperatures that are 2-3 orders of ma gnit ude 
larger than in simple metals and most alloys ]2(i gjj. 

In previous neutron scattering experiments fast dy- 
namics in viscous Zr46.8Tig.2Cu7.5NiioBe27.5 have been 
investigated: Data analysis in the framework of MCT is 
in full agreement with the predicted scaling functions and 
gives a T c some 20 % below Tu q . In an experiment on 
viscous Pd4oNiioCu3oP20 (Pd40) the focus was on slow 
dynamics [p3| : The long-time decay of correlations ex- 
hibits the common features of glass forming liquids, i.e. 
structural relaxation obeys a stretching in time and a 
universal time-temperature superposition. 

Here, inelastic neutron scattering results on liquid 
Pd4 3 Ni 10 Cu 2 7P2o (Pd43) are reported. The small change 
in composition compared to the Pd4 Ni 10 Cu3oP2o results 
in an improved stability with respect to crystallization 
|l6f - a cooling rate as low as 0.09 K/s is sufficient 
to avoid crystallization and to form bulk metallic glass 
This allows measurements of transport coefficients 
at temperatures around the mode coupling T c with ex- 
perimental techniques that only require heat treatment 
for some minutes, e.g. rheometry [p0[ or radio tracer dif- 
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Compared to previous neutron scattering measure- 
ments |22| , p3| a different experimental setup has been 
used that covers an extended range of momentum trans- 
fers q. This allows a detailed investigation of the fast 
MCT j3 relaxation and the correlation between struc- 
tural relaxation and long range atomic transport. In 
addition, the temperature range has been extended by 
several 100 K up to temperatures at which the diffusion 
of the atoms is of the order of that in simple monoatomic 
metals. The results are set in context to macroscopic 
transport coefficients, the dynamics in simple liquids and 
the predictions by the mode coupling theory of the liquid 
to glass transition. 



II. EXPERIMENTAL 

Pd43NiioCu27P2o w & s prepared from a mixture of pure 
elements by induction melting in a silica tube. The melt 
was subject to a B2O3 flux treatment in order to remove 
oxide impurities. Differential scanning calorimetry with 
a heating rate of 40 K/min resulted in a T g at 578 K and a 
Tu q at 863 K in accordance with literature values |l^, ^4| . 
For the neutron time-of-ffight experiment a thin-walled 
AI2O3 container has been used giving a hollow cylinder 
sample geometry of 22 mm in diameter and a thickness of 
1 mm. During the measurement the liquid was covered by 
a thin layer of n B203 flux material. For the chosen ge- 
ometry and neutron wavelength the sample scatters less 
than 2%. Therefore, multiple scattering, which would 
alter the data especially towards low q |2(|, only has a 
negligible effect. 

Microscopic dynamics in liquid PdNiCuP has been 
investigated on the neutron time-of-ffight spectrometer 
IN 6 at the Institut Laue-Langevin in Grenoble. An 



incident neutron wavelength of A = 5.1 A yielded an 
accessible wavenumber range at zero energy transfer of 
q = 0.75 — 1.95 A at an energy resolution of 92 fj,eV 
(FWHM). Regarding the scattering cross sections of the 
individual elements PdNiCuP is a 90% coherent scat- 
terer. However, with the first structure factor maximum 
at go — 2. 9 A -1 our spectra are dominated by incoherent 
scattering, that is dominated by the contributions from 
Ni with ~ 73 % and Cu with ~ 20 %. 

The alloy was measured at room temperature to obtain 
the instrumental energy resolution function. Data were 
collected in the liquid between 833 K and 913 K in steps of 
40 K and between 973 K and 1373 K in steps of 100 K with 
a duration between 2 and 5 hours each. At each tempera- 
ture empty cell runs were performed. The data at 833 K, 
30 K below the liquidus, do not show signs of crystalliza- 
tion, which is in accordance with the time-temperature 
transformation of the undercooled liquid ]24j ] . In order to 
obtain the scattering law S(q, u) (Fig. |l|), raw data were 
normalized to a Vanadium standard, corrected for self- 
absorption and container scattering, interpolated to con- 
stant q, and symmetrized with respect to energy with the 
detailed balance factor. Fourier deconvolution of S(q, u>) 
and normalization to 1 for t = gives the self correlation 
function t). 



III. MODE-COUPLING THEORY 
ASYMPTOTICS 

The mode-coupling theory of the liquid to glass tran- 
sition (MCT) 1 11, [lj] is developed in the well-defined 
frame of molecular hydrodynamics p7j . Starting with 
the particle density, p r (t) — i- e - the po- 

sitions, Ri, of each particle, i, at time t, the Zwanzig-Mori 
formalism provides an exact equation of motion for the 
density correlation function, & q (t) = (p q {t)*p q ) / (|p 9 | 2 ): 



*,(*) + n» $,(*)+ / M q (t-t')$ g (t')dt' = (i) 

Jo 

where fl q represents a phonon frequency at wavenumber 
q and M q (t) a correlation function of force fluctuations 
which in turn is a functional of density correlations. The 
static part of the fluctuating force is a linear combination 
of density fluctuation pairs and depends therefore only on 
the interaction potential. 

Since MCT does not aim to describe the detailed mi- 
croscopic phonon distribution, but rather to give a uni- 
versal picture of relaxational dynamics at longer times, 
the integral kernel is split into 



M q (t) = v q 6(t) +{l 2 q m g (t), 



(2) 



where v q models the damping by "fast" modes and m q (t) 
accounts for memory effects through the coupling of 
"slow" modes. The basic idea of the mode-coupling the- 
ory of the liquid to glass transition is to consider as 
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"slow" all products of density fluctuations. By deriva- 
tion, m q (t) contains no terms linear in $ q (i). Therefore, 
in lowest order, it is a quadratic functional 



m q (t) 



E 

<?l+<?2 = 



V q ( qi ,q 2 )^ qi (t) $ q2 (t). 



(3) 



In this approximation, the coupling coefficients V q are 
specified in terms of the static structure of the liquid. 

Eqs. 1, 2 and 3 lead to the following scenario: a fast (3- 
relaxation process, which can be visualized as a rattling 
of the atoms in the cages formed by their neighbouring 
atoms, prepares structural a relaxation, responsible for 
viscous flow. At an ideal glass transition temperature T c 
the transport mechanism crosses over from glass-like ac- 
tivated hopping processes to liquid-like collective motion. 
In other words, at T c , the cages are no longer stable on 
the time scale of a diffusive jump. 

A common feature of structural a relaxation in glass- 
forming liquids is a stretching in correlation functions 
over a wider time range than expected for exponential 
relaxation |], [TJJ] . Experimental data in the a relaxation 
regime can usually be well described by a stretched ex- 
ponential function 



F(q,t) = f c q exp[-(t/r q f"} 



(4) 



with an exponent [3 q < 1. r q is the relaxation time and 
f q < 1 accounts for the initial decay of correlations due 
to phonons and the fast relaxation process. For tempera- 
tures above T c mode coupling theory predicts a universal 
time-temperature superpostion of structural relaxation 
resulting in a temperature independent stretching expo- 
nent j3 q . Many glass- forming liquids in contrast exhibit a 
temperature dependence of (3 — in some (3 is decreasing, 
in others increasing on temperature increase. 

Solutions of the mode coupling equations for a hard 
sphere system confirm, that the MCT a relaxation is 
well described by a stretched exponential function and 
Eq. (H) becomes a special, long time solution of mode 
coupling theory p3[ . For the mean relaxation times 
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(r q ) = J dtF(q,t)/f° = Tq pr x T{pr x ) 
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mode coupling theory predicts that (r q {T)) is inversely 
proportional to the diffusivity D(T) and obeys an asymp- 
totic scaling law for temperatures above but close to T c : 



(r q ) cx 1/D cx [(T-r c )/T c ; 



(6) 



Asymptotic expansions of Eq. (0^) show that in the in- 
termediate /3-relaxation regime around a crossover time, 
t a , and for temperatures close to T c , the asymptotic form 
of the correlation function is independent of the detailed 
structure of the coupling coefficients and exhibits a uni- 
versal factorization property: 
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FIG. 1: Scattering law S(q,u>) (logarithmic scale) of liq- 
uid Pd43NiioCu27P20 obtained on the neutron time-of-flight 
spectrometer IN 6. The signal is dominated by the incoherent 
scattering from Ni and Cu. The data at 300 K represent the 
instrumental energy resolution. Measurements in the liquid 
(Tu q = 863 K) reveals a broad quasielastic signal. 



where f q represents the Debye- Waller factor and h q an 
amplitude. The scaling function g\(t) is defined by just 
one shape parameter A. Close to T c mode coupling theory 
predicts the temperature dependence of t a and h q with 
the asymptotic scaling functions: 

t a oc (T - T c y 1/2a and h q cx (T - T c ) 1/2 . (8) 

There are tables providing a, 7 and g\ as a function of 
A M. 



IV. RESULTS 

Figure [I] shows the scattering law S(q,oj) of liquid 
PdNiCuP at q = 1.05 A -1 . The signal is dominated by the 
incoherent scattering from Ni and Cu. S(q, u>) displays a 
quasielastic line with an increasing width on temperature 
increase and with wings extending up to several meV. 
Above some 3 meV the quasielastic signal merges into a 
constant as expected for an ideal Debye solid and found 
in liquid GeC>2 [|9|. The quasielastic signal of S(q,u>) is 
fairly well separated from vibrations. 



A. Vibrations 

There is no general understanding of the influence of 
phonons on the atomic transport in liquids. Whereas 
in monoatomic alkali melts one finds low-lyingcollective 
phonon modes that mediate mass transport ||, [30[ [H], 
the mode coupling theory of the liquid to glass transition 
does not consider the detailed phonon dynamics for its 
description of atomic transport in viscous liquids (Eq. 0). 
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FIG. 2: Density of states g(uj) representing mainly the Ni and 
Cu vibrations in liquid Pd43NiioCu27P20- With a maximum 
at some 17meV (corresponding to some 0.04 ps) vibrations 
are fairly well separated from the fast relaxational dynamics 
around 1 ps. The small expansion coefficient reflects itself in 
a weak temperature dependence of g(oj). 



In particular, dynamics in hard sphere- like colloidal sus- 
pensions, that do not exhibit vibrations, are in excellent 
agreement with MCT predictions |J. 

In most molecular liquids broad phonon distributions 
overlap with the quasielastic signal. The phonon density 
of states exhibits a first maximum usually between 4 and 
8meV H|, |33| . This aggravates the theoretical descrip- 
tion of the relaxational dynamics. Figure ^| displays the 
vibrational density of states g{oS) representing mainly the 
incoherent contributions of the Ni and Cu atoms of liq- 
uid PdNiCuP. g(ui) has been derived from the scattering 
law S(q,u>) using a procedure that assumes pure inco- 
herent scatte ring and uses the multiphonon correction as 
described in |34| . 

g{ui) displays a weak temperature dependence indi- 
cating that anharmonic contributions to the interatomic 
potential are minor. This is in accordance with an ex- 
pansion coefficient of only 4 x 10 _5 K _1 p9| , that is an 
order of magnitude smaller than that of most molecu- 
lar glass forming liquids. With a maximum in g(u>) at 
some 17meV, vibrations are fairly well separated from 
the quasielastic signal that extends up to several meV. 
In addition, PdNiCuP does not exhibit a "boson peak" 
- a maximum in S(q, ui) that is found in other glasses 
usually at a few meV - in the q range investigated. This 
allows a detailed investigation of the atomic transport in 
liquid PdNiCuP. 



obtained by Fourier transformation of measured S(q,u>), 
division of the instrumental resolution function, and nor- 
malization with the value at t = 0. Between and ~ 1 ps 
phonons and a fast process lead to a decrease in t) 
from 1 towards a plateau. Figure [| displays the long-time 
decay of $(<7, t) at q= 1.05 A -1 from this plateau towards 
zero in a semilogarithmic representation. 

In contrast to simple metallic liquids, liquid PdNiCuP 
exhibits a structural relaxation that shows stretching in 
time. The lines are fits with the stretched exponential 
function (Eq. In a first fitting process the stretching 
exponent (3 q (T) was treated as a parameter: j3 q (T) varies 
weakly around a mean (3 — 0.75. For the further data 
analysis a /3 = 0.75 was used. 

Structural relaxation described by mode coupling the- 
ory displays stretching and a time-temperature super- 
position. Using results from the fitting procedure with 
Eq. ||and a stretching exponent /3 = 0.75 master curves 
$>(q,(t/T q ))/f q have been constructed for the data in 
the structural relaxation regime above some 1 ps. Fig- 
ure U shows rescaled (t/T q ))/f q at q = 1.05 A -1 . 
Over the entire temperature range, that spans 540 K, 
the <§>{q,{t/T q ))/ f q fall on a master curve: a time- 
temperature superposition of structural relaxation holds. 
We note, that the validity of the time-temperature su- 
perposition up to temperatures at which th e tran sport 
coefficients approach that of simple liquids (IVC) is in 
marked contrast to the idea, that viscous liquids exhibit 
a transition to a liquid with the stretching exponent (3 
approaching 1 on temperature increase. 

Stretching of self correlation functions is generally 
found to be more pronounced in fragile glass-forming liq- 
uids, characterized by a curved temperature dependence 
of viscosity in an Arrhenius plot |3q]; e.g. the van der 
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B. Structural Relaxation 

The quasielastic signal is best analyzed in the time do- 
main with a removal of the instrumental resolution func- 
tion. The density correlation function <I>(q,i) has been 



FIG. 3: Normalized density correlation function <fr(q,t) of 
liquid PdNiCuP at 1.05 A -1 as obtained by Fourier decon- 
volution of measured S(q,ui). Structural relaxation causes 
the final decay of &(q,t) — > 0. The lines are fits a stretched 
exponential function (Eq. ^). 
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FIG. 4: Rescaling of the density correlation function in the 
a relaxation regime (for t > 1 ps) using results from fits with 
a stretched exponential function: a time-temperature super- 
position of structural relaxation holds from Tu q — 30 K to 
Tuq + 510 K. The line is a fit with Eq. ^| resulting in a stretch- 
ing exponent p — 0.75. 
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FIG. 5: Rescaling of the density correlation function $(g, t) 
at 1073 K. Within the accessible q range the stretching of 
the structural relaxation is fairly q independent. The line 
represents a stretched exponential function (Eq. ^) with a 
stretching exponent /3 = 0.75. 



Waals liquid orthoterphenyl with a /3~0.5 pjj\ . In an in- 
termediate system like hydrogen-bond forming glycerol 
a (3 ~ 0.6 has been reported |56). The stretching ex- 
ponent in liquid PdNiCuP metals compares well to the 
value /3~0.75 found in covalent-network forming sodium 
disilicate melts |37j and indicates that in this context 
the PdNiCuP alloy might classify as a fairly strong glass 
forming liquid. 

Figure | shows rescaled $(q,(t/T q ))/f q at 1073 K for 
q values in the range between 0.75A -1 and 1.95A -1 . 
Structural relaxation in liquid PdNiCuP can be described 



with a stretched exponential function and a q and tem- 
perature independent stretching exponent (3 — 0.75 ± 
0.02. [3 q shows a small but systematic decrease with 
increasing q in agreement with MCT calculations for a 
hard-sphere system j^]. However, a small variation of 
f3 q in the fitting procedure has no significant effect on 
the resulting mean relaxation times (Eq. |^). 



C. Diffusion 

In the hydrodynamic limit for q — > one expects that 
the mean relaxation times (t 9 ) are indirect proportional 
to the square of the momentum transfer q [p7[ . In liq- 
uid PdNiCuP the 1 / q 2 dependence holds even up to 
1.9 A -1 (Fig. ||) while the first structure factor maxi- 
mum is at qo ~ 2.9 A -1 . In liquid alkali-metals one ob- 
serves a systematic deviation from a q 2 dependence of 
the quasielastic line width already at q values that cor- 
respond to 1/10 of their structure factor maximum qo- 
This deviation is explained by low-lying phonon modes 
that mediate atomic transport || |3l] ]. Mean relaxation 
times for self motion in the MCT hard sphere system, 
in contrast, vary rather well proportionally to 1 / q 2 for 
q values extending even above the first structure factor 
maximum p"3] ] . We note that mixing of hard spheres with 
different radii should enforce this behaviour even more. 
It appears that the validity of 1 / (r q ) oc q 2 also for inter- 
mediate q values is a signature of the atomic transport 
mechanism in dense liquids. 

The 1 / q 2 dependence of the mean relaxation times, 
also demonstrates that structural relaxation leads to long 
range atomic transport with a diffusivity f27j] 
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FIG. 6: Mean relaxation times (T q ) of tagged particle motion 
in liquid Pd43NiioCu27P20- 1 / {T q } shows a q 2 dependence as 
expected for long range atomic transport for q — » 0. The slope 
corresponds to the self diffusion coefficient D(T). 
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FIG. 7: Self diffusion coefficient D of Ni and Cu in liquid 
Pd43NiioCu27P20 (closed circles) and Pd4oNiioCu3oP20 (open 
circles from Ref. derived from the mean relaxation times. 
The lines are fits with MCT r scaling (Eq. ^|). The data are 
well described with T c = 700 ± 30 and 7 = 2.7 ± 0.2. 
For c omp arison 57 Co tracer diffusion in Pd43NiioCu27P20 
(Ref. Ill]) and diffusivities calculated from viscosity T) of liq- 
uid Pd4oNiinCu3oP20 (Ref. ^J) via the Stokes-Einstein rela- 
tion (Eq. are shown. 



Figure displays the diffusivities D in liquid Pd43Nii - 
CU27P20 as a function of 1 / T. Values range from 
3 ± lx 10 _n m 2 s _1 at 833K to 1.4 ± 0.3 x lO^nvV 1 
at 1373 K. At the liquidus diffusion is about 2 orders of 
magnitude slower as compared to simple metallic liquids 
and most alloys. Diffusivities in liquid Pd40 show a simi- 
lar temperature dependence but are larger by some 20 %. 
The smaller mobility of the atoms in Pd43 appears to 
come along with the slightly better glass forming ability 
of the Pd43 alloy. 

The lines in Figure fj] represent fits with the MCT r 
scaling law (Eq. ||) to the diffusivities D(T). Although 
Eq. H is only valid close to T c , Eq. || allows a rough esti- 
mate of the crossover temperature T c and the exponent 
7. A fit to D(T) of liquid Pd43 at temperatures up to 
1073K yields a T C = 700±30K and a 7 = 2.7±0.2. The 
exponent 7 compares well to the 7 = 2.62 found in the 
MCT hard sphere system [jlj|. The temperature depen- 
dence of the diffusivity in liquid PdNiCuP alloys is dif- 
ferent from the D oc T n behaviour (with n ~ 2) expected 
for uncorrelated binary collisions of hard spheres Jl], Q 
and found with inelastic neutron scattering in expanded 
liquid alkali- metals at high temperature |39| . 

Convection effects are a severe problem in macroscopic 
diffusion measurements in ordinary liquids under gravity 
conditions with respect to the absolute value of the self 
diffusivity D and its temperature dependence. A fig ex- 
periment on liquid Sn revealed that well above the melt- 
ing point convection even dominates the mass transport 
p0[ . The fig data in liquid Sn obey D oc T 2 for tem- 
peratures between T m and 2 x T m . Liquid Sn exhibits a 



packing fraction of ~ 0.4 as in liquid alkali-metals. Inelas- 
tic neutron scattering data are not affected by convection 
because it probes dynamics on significantly shorter times. 

Diffusivities from recent 57 Co tracer diffusion measure- 
ments in liquid Pd43 p5| ] are in excellent agreement with 
the diffusivities obtained from inelastic neutron scatter- 
ing (Fig. ^). This indicates that the mobility of Co is 
very similar to the Ni and Cu mobility observed in the 
incoherent neutron scattering signal and that convection 
effects play no significant role in the tracer diffusion ex- 
periment. The latter is in accordance with a viscosity 77 
that is two orders of magnitude larger as compared to 
simple metallic melts J2l| . 

In simple liquids shear viscosity rj and the self dif- 
fusivity D generally obey the Stokes-Einstein relation 

Sill: 



D = kbT/ {(Sirarf), 



(10) 



with reasonable values for the hydrodynamic radius a. 
Eq. [h] even holds in most molecular liquids jllj . Figure 
shows the diffusivity D v calculated from the viscos- 
ity data of liquid Pd40 §| via Eq. |l0] using a = 1.15 A 
that represents the Ni hard sphere radius (Cu: 1.17 A) 
from Ref. The Stokes-Einstein relation holds with 
D v = D. We note, that a is similar to the mean next 
nearest neighbor distance displayed in the static struc- 
ture factor d = 2n/qo ~ 1.1 A. In multicomponent liq- 
uids, that exhibit a dense packing of hard spheres with 
comparable hard sphere radii one expects a similar mo- 
bility of the different components. In PdNiCuP the hard 
sphere radii of the different atoms have values within 
20 %. Because viscous flow represents the dynamics of 
all components, the validity of Eq. [l(] in liquid PdNiCuP 
demonstrates that above Tu q the mobility of the large 
and numerous Pd atoms is quite similar to the Ni and 
Cu atoms ||. 

In colloidal suspensions the mean relaxation times 
(r)((p) and the inverse of the diffusivity 1/D(ip) show the 
same dependence as a function of the packing fraction. 
The slope for (r((/?)) and 1/D((p) results in a 7 = 2.7 

f. In molecular dynamics simulations on viscous NiP 
in contrast, the asymptotic MCT t scaling prediction 
(Eq. |) is violated: (r(T)) and 1/D(T) do not exhibit the 
same temperature dependence. Figures and [l^ demon- 
strate that in liquid PdNiCuP Eq. || holds quite well. 



D. Fast (3 relaxation 

Within the mode coupling theory of the liquid to glass 
transition, structural relaxation and therefore long range 
mass transport is preceded by a fast (3 relaxation process 
whose time scale is typically in the order of a picosecond. 
For q values well below the structure factor maximum 
and incoherent scattering the amplitude of the (3 relax- 
ation is increasing with increasing q. Figure ^| displays 
the density correlation function of liquid PdNiCuP at 
q = 1.95 A -1 . Between and ~ 0.2 ps phonons lead to 
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a rapid decay of atomic correlations (represented by the 
first data point in &(q,t) in Fig. ||), so that t) de- 
creases from 1 towards a plateau. On approaching this 
plateau a fast relaxation is clearly seen in <5>(q, t) through 
an additional intensity below some 1 ps. The lines are fits 
with the MCT (3 scaling law (Eq. @). For q < 1.4 A" 1 the 
limited dynamic range of the instrument prevents data 
analysis in the fast j3 relaxation regime. Above 973 K 
$(</, t) can not consistently be described with the asymp- 
totic scaling laws (Eq. |t],|§)- The dynamic range in which 
Eq. |?] holds is increasing with decreasing temperature. 

The data were fitted in a two-step procedure: starting 
with an arbitrary line shape parameter A, fits to individ- 
ual curves were used to estimate the scaling factors f q , 
h q , and t a which physically represent the Debye- Waller 
factor, the amplitude of /3-relaxation, and a characteris- 
tic time of /3-relaxation. Using these values, the &(q,t) 
measured at different temperatures were superimposed 
on to master curves ($>(q,t/t a ) — f q )/h q . After fixing a 
q- independent mean t a , the fit yielded a temperature and 
q independent A = 0.78 ± 0.04. 

This result compares well to the hard spheres value 
A = 0.766 of the numerical MCT solution ||], to the 
A = 0.77 ±0.04 found in liquid Zr46.8Tig.2Cu7.5NiioBe27.5 
[22 1 and it is similar to the values in other viscous liquids 
[11, B3. A = 0.78 defines an exponent of the r scaling 
law (Eq. ||) of 7 = 2.7 psfl . This is consistent with the 
temperature dependence of the mean relaxation times 
that give 7 = 2.7 ±0.2 (Fig. @). 

Figure [)] shows amplitude h q and time scale t a of the 
fast j3 relaxation in liquid PdNiCuP rectified according 
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FIG. 9: Mean amplitude (h q ) (triangles) and time scale t a 
(circles) of fast j3 relaxation in liquid Pd43NiioCu27P20 rec- 
tified according the MCT predictions (Eq. ^), with an expo- 
nent a = 0.29 defined by A = 0.78. The inset shows the q 
dependence of the amplitude h q between 833 K (•) and 953 K 
(o). Both h q and t a extrapolate to T c = 720±20K. A and 
T c obtained from the analysis of the (3 relaxation regime are 
consistent with T c = 700±30K and 7 = 2.7±0.2 from the 
analysis of the mean a relaxation times (Eq. ^). 



is in accordance with the MCT predictions. Both h q 
and t a extrapolate to T c = 720±20K that is close to 
the T c = 700±30K obtained via the r scaling law (Eq. 
||). The relaxational dynamics in liquid Pd43Nii Cu27P20 
can consistently be described with the universal scaling 
functions of the mode coupling theory of the liquid to 
glass transition at temperatures up to ~ 1.3 x T c with 
T C = 710K and A = 0.78. 
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FIG. 8: Normalized time correlation function $(g, t) of liq- 
uid Pd43NiioCu27P20 in the /3 relaxation regime around one 
picosecond. The solid lines represent fits with the MCT /3 
scaling function (Eq. [?]) giving a temperature and q indepen- 
dent shape parameter A = 0.78 ±0.04. Above ~ lps data 
are well described by a stretched exponential function (Eq. ^| 
dashed lines). 



E. Atomic transport mechanism 

Pd43Nii Cu27P2o exhibits an excellent glass forming 
ability: a cooling rate of only 0.09 K/s is sufficient to 
avoid crystallization and to form a glass |24f| . Conse- 
quently transport coefficients can continuously be mea- 
sured from the equilibrium liquid down to the conven- 
tional glass transition temperature T g . Figure |lC| displays 
diffusivitics in viscous PdNiCuP alloys from Co tracer 
diffusion measurements, that cover more than 13 decades 
p5| . p4| , derived from the mean relaxation times and from 
viscosity jl9[ |l) calculated via Stokes-Einstein (Eq. [To]) . 

At 1373 K the diffusivity in liquid PdNiCuP is 1.4 x 
10 _9 m 2 s _1 which is similar to diffusivities known from 
monoatomic liquid metals at their respective melting 
temperature. E.g. close to the melting point inelastic 
neutron scattering reveals D — 4.2 x 10 _9 m 2 s _1 in 
liquid sodium Q and tracer diffusion under fj,g condi- 
tions D = 2.0 x 10- 9 m 2 s" 1 in liquid tin @. At high 
temperatures the Stokes-Einstein relation holds in liquid 
PdNiCuP (Fig. Q). On lowering the temperature diffu- 
sion slows down drastically. Time scales for viscous flow 
and Co tracer diffusion start to decouple in the vicinity 
of the mode coupling T c ~ 710 K and differ already by 
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more than 3 orders of magnitude at 600 K. Also around T c 
Co tracer diffusivities merge into an Arrhenius-type tem- 
perature dependence (D(T) = D exp {—H/UbT), where 
Dq is a prefactor and H the activation enthalpy) that 
extends down to the glass transition |25|, 

The temperature dependence of the transport coef- 
ficients in Figure strongly supports the MCT pre- 
diction of a change in the atomic transport mechanism 
from viscous flow to glass like hopping at T c . There 
are extensive tracer diffusion measurements in super- 
cooled Zr46.8Tig.2Cu7.5NiioBe27.5 at temperatures be- 
tween the glass transition and about 200 K below T c 
po| . The diffusivities of various tracers exhibit a size de- 
pendence: the smaller the atoms the faster they diffuse 
and the smaller the activation enthalpy H. In viscous 
Zr46.8Tig.2Cu7.5NiioBe27.5 the diffusivities of the various 
atoms approach each other with increasing temperature. 
Because in this alloy crystallization prevents access to the 
temperature range around T c , one can not conclude from 
the present data whether the diffusivities of the various 
tracers merge at the MCT T c . 

The isotope effect E = (D a /Df, — l)/(\/mb/m a — I) 
(diffusivity D and mass m of two isotopes a, b) is a mea- 
sure of the degree of collectivity of the atomic transport. 
For diffusion via single jumps in dense packed lattices E 
is generally in the order of unity |Q. For uncorrelated 
binary collisions one also expects E — > 1 . A vanishing iso- 
tope effect indicates a collective transport mechanism in- 
volving a large number of atoms. In liquid Sn an isotope 
effect of about 0.4 has been reported, that is increasing 
with increasing temperature |40|. MD simulations on a 
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FIG. 10: Mass transport in PdNiCuP melts: At high tem- 
peratures diffusivities are similar to diffusivities in liquid Sn 
and Na close to T m (marked by the arrows). Time scales for 
viscous flow and Co tracer diffusion start to decouple and Co 
tracer diffusivities merge into an Arrhenius-type temperature 
dependence that extends down to the glass transition. Both 
occurs in the vicinity of mode coupling T c . The line represents 
the MCT t scaling law with T C = 710K and 7 = 2.7. 



binary Lennard-Jones liquid demonstrate that changes in 
the density in the order of 20 % only result in a continu- 
ous increase in E from to ~ 0.3 [47|] . In metallic glasses 
[ 48 1 as well as in supercooled Zr4 6 . 8 Ti 8 .2Cu7.5Nii Be27.5 
[ 49 j and Pd4oNiioCu3oP2o |4^] the isotope effect is close 
to zero as a result of a highly collective, thermally ac- 
tivated hopping process. Measurements of the isotope 
effect in Pd 43 Nii Cu27P20 H give E ~ 0.05 over the 
entire temperature range from the glass transition to the 
equilibrium liquid. 

The atomic transport mechanism in dense PdNiCuP 
remains highly collective even in the equilibrium liquid 
in contrast to the findings in simple metallic liquids. 
Viscous PdNiCuP exhibits an expansion coefficient of 
4 x lO^K" 1 @. Between 1373K and 833K, the den- 
sity and therefore, assuming a temperature independent 
hard sphere radius, the packing fraction changes only by 
about 2 %. An increase in the radius of the attributed 
hard sphere radii Q with increasing temperature results 
in an even smaller change in the packing fraction. This 
small increase in the packing fraction causes the drasti- 
cal slowing down of dynamics on temperature decrease 
in accordance with the mode coupling scenario. Dense 
packing and the resulting collective transport mechanism 
extends also to high temperatures, where diffusivities are 
similar to that in simple metallic liquids. This indicates 
that the atomic transport mechanism in liquids is not 
controled by the value of the transport coefficients but 
rather by that of the packing fraction. 



V. CONCLUSIONS 

Atomic transport in liquid Pd43NiioCu27P2o bas been 
investigated with inelastic, incoherent neutron scatter- 
ing. The PdNiCuP melt is characterized by a packing 
fraction of about 0.52 that is some 20 % larger as com- 
pared to monoatomic alkali-melts. The self correlation 
function shows a two-step decay as known from other 
non-metallic glass- forming liquids: a fast relaxation pro- 
cess precedes structural relaxation. The structural re- 
laxation exhibits stretching with a stretching exponent 
(3 ~ 0.75 and a time-temperature superpostion that holds 
for temperatures as high as 1.75 x T c . The relaxational 
dynamics can consistently be described within the frame- 
work of the mode coupling theory of the liquid to glass 
transition with a temperature and q independent line 
shape parameter A ~ 0.78. Universal MCT scaling laws 
extrapolate to a critical temperature T c ~ 710 K some 
20% below the liquidus. 

Diffusivities derived from the mean relaxation times 
compare well with Co diffusivities from tracer diffusion 
measurements. Above T c diffusivities calculated from vis- 
cosity [^l| via the Stokes-Einstein relation are in excel- 
lent agreement with the diffusivities measred by tracer 
diffusion [ p5[ and neutron scattering. In contrast to sim- 
ple metallic liquids the atomic transport in dense liquid 
PdNiCuP is characterized by a drastical slowing down of 
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dynamics on approaching T c , a q~ 2 dependence of the 
mean relaxation times at intermediate q and a vanishing 
isotope effect as a result of a highly collective transport 
mechanism in the dense packed liquid. At temperatures 
as high as 2xT c diffusion in liquid PdNiCuP is as fast 
as in simple monoatomic liquids at their melting points. 
However, the difference in the underlying atomic trans- 
port mechanism indicates that the diffusion mechanism 
in liquids is not controlled by the value of the diffusivity 
but rather by that of the packing fraction. 
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